3530  IEEE  TRANSACTIONS  ON  ANTENNAS  AND  PROPAGATION,  VOL.  57,  NO.  1 1,  NOVEMBER  2009 

Estimating  Moving  Targets  Behind 
Reinforced  Walls  Using  Radar 

Marija  M.  Nikolic.  Arye  Nehorai,  Fellow,  IEEE ,  and  Antonije  R.  Djordjevic 


Abstract — We  consider  the  estimation  of  moving  targets  located 
behind  concrete  walls  reinforced  with  metallic  bars,  using  radar 
measurements.  The  periodic  structure  of  the  rebar  severely  atten¬ 
uates  and  distorts  transmitted  waveforms,  which  produces  defo- 
cused  images  with  ghost  target  estimates.  We  apply  beamforming 
to  estimate  permittivity  and  thickness  of  the  wall  and  number  and 
position  of  the  targets.  The  proposed  solution  is  based  on  accurate 
physical  models  calculated  using  the  method  of  moments.  We  show 
that  the  estimation  is  significantly  improved  by  modeling  the  wave¬ 
form  distortion  due  to  the  bars.  The  resulting  images  are  focused 
and  clearly  represent  the  contours  of  the  targets.  The  algorithm  is 
robust  to  the  ambiguities  in  bar  parameter  values.  In  addition,  the 
minimal  necessary  SNRis  lower  compared  with  the  case  in  which 
the  influence  of  the  bars  on  the  signal  shape  is  ignored. 

Index  Terms — Beamforming,  periodic  structure,  reinforced  wall, 
through-the-wall  sensing. 


I.  Introduction 

AN  important  task  in  urban  warfare  is  using  exterior 
electromagnetic  sensing  to  detect  people  hidden  inside 
buildings.  Solving  this  problem  is  challenging  because  of  the 
complex  and  unknown  environment.  The  walls  significantly 
attenuate  and  distort  transmitted  signals;  hence  electromagnetic 
modeling  has  an  important  influence  on  estimation  accuracy. 
The  distortion  is  particularly  prominent  in  the  case  of  walls 
reinforced  by  parallel  steel  bars  (rods)  or  square-grid  meshes, 
which  are  commonly  used  in  construction  [8]— [  1 6] .  We  show 
that  the  complex,  periodic  structure  of  reinforced  walls  imposes 
significant  difficulties  in  the  target  imaging  that  are  not  encoun¬ 
tered  in  the  case  of  homogenous  walls.  If  the  influence  of  the 
bars  is  not  taken  into  account,  images  of  the  objects  located 
behind  the  wall  become  defocused  and  false  targets  appear. 

Several  studies  considered  through-the-wall  estimation 
[1]— [7],  In  [1]  and  [2]  the  wall  is  treated  as  a  homogeneous 
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dielectric  slab  of  known  thickness  and  permittivity.  The  case  in 
which  the  wall  parameters  are  unknown  is  studied  in  [3].  The 
authors  of  [3]  exploit  measurements  from  two  or  more  standoff 
distances  from  the  wall  to  calculate  images  for  various  wall 
parameters.  Assuming  that  the  most  focused  image  is  obtained 
using  accurate  wall  parameters,  an  image-focusing  metric  is 
applied  to  estimate  the  correct  wall  thickness  and  dielectric 
permittivity.  Imaging  through  walls  made  of  cinder  block  is 
addressed  in  [6].  The  authors  of  this  article  consider  the  wall  as 
a  periodic  structure  and  apply  the  SAR  algorithm  to  estimate  the 
target  position.  However,  besides  true  target  position,  the  SAR 
algorithm  also  produces  a  ghost  target  image.  An  estimation  of 
the  building  layout  is  conducted  in  [7]. 

Electromagnetic  propagation  through  reinforced  concrete 
walls  is  well  covered  in  the  literature  for  the  purpose  of  mobile 
channel  modeling.  Rigorous  analysis  of  the  transmission  and 
reflection  coefficients  of  reinforced  walls  has  been  carried  out, 
e.g.,  in  [8] — [10].  A  non-destructive  analysis  of  reinforced  con¬ 
crete  using  electromagnetic  waves  was  conducted  in  [1 1] — [14]. 
Here,  our  focus  is  on  estimating  the  position  of  the  targets 
located  behind  reinforced  walls.  We  separately  investigate 
the  cases  of  known  and  unknown  wall  parameters.  We  apply 
beamforming  to  the  wideband  measurements  made  by  a  sensor 
array  to  estimate  wall  thickness  and  dielectric  permittivity.  The 
number  and  positions  of  the  targets  are  further  estimated  using 
beamforming  [1],  as  well  as  wall-parameter  estimates.  In  the 
case  in  which  signal  distortion  due  to  the  bars  is  not  taken  into 
account,  the  target  spreads  are  large.  We  show  that  significant 
improvement  is  achieved  when  bar  characteristics  are  partly 
or  completely  known.  This  improvement  is  particularly  pro¬ 
nounced  when  the  signal-to-noise  ratio  is  low. 

The  organization  of  the  paper  is  as  follows.  In  Section  II,  we 
discuss  the  signal  propagation  through  reinforced  walls.  We  also 
present  measurement  models  and  our  electromagnetic  solver. 
Estimation  of  wall  permittivity  and  thickness  is  discussed  in 
Section  III.  In  Section  IV,  we  study  moving-target  estimation 
with  unknown  and  known  wall  parameters.  Finally,  in  Section  V 
we  present  simulation  results. 

II.  Signal  Propagation  Through  Reinforced  Walls 
A.  Reinforced  Walls 

Reinforced  walls  are  generally  one  of  two  types:  with  bars  or 
with  mesh.  When  the  incident  electromagnetic  field  is  parallel 
to  the  vertical  bars,  the  horizontal  bars  (cross  bars)  have  small 
influence  on  the  transmission  and  reflection  coefficient  of  the 
wall  [8].  The  same  reasoning  applies  when  the  incident  electric 
field  is  parallel  to  the  horizontal  bars.  Therefore,  without  loss 
of  generality,  we  consider  a  wall  reinforced  with  vertical  bars. 
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Fig.  1 .  Cross  section  of  reinforced  concrete  wall. 


We  depict  the  cross  section  of  such  a  wall  in  Fig.  1,  where 
is  the  distance  between  adjacent  bars  (bar  period),  D),;lr  is  the 
bar  diameter,  and  w  is  the  wall  thickness.  We  assume  the  bars  to 
be  perfectly  conducting  and  immersed  in  a  homogeneous  ma¬ 
terial  (concrete)  of  relative  dielectric  permittivity  cr.  We  con¬ 
sider  a  two-dimensional  problem,  since  electromagnetic  mod¬ 
eling  of  electrically  large  structures,  such  as  walls,  is  extremely 
time  consuming. 
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Fig.  2.  Signals  reflected  from  a  perfectly  conducting  plane  (PEC)  located  be¬ 
hind  a  reinforced  wall  and  a  homogeneous  wall. 


B.  Measurement  Model 

We  determine  the  frequency  response  of  the  considered 
system  using  a  surface  formulation  with  the  method  of  mo¬ 
ments  [17],  [18],  The  program  that  we  developed  uses  the 
equivalence  principle  to  divide  the  system  under  consideration 
into  a  number  of  subsystems  (entities),  each  filled  with  a  ho¬ 
mogeneous  medium.  In  our  case,  the  system  consists  of  sensors 
(line  conductors  whose  cross  section  is  electrically  small)  and 
dielectric  walls  with  metallic  bars  inside.  The  program  excites 
one  sensor  at  a  time  by  an  impressed  electric  field,  calculates 
the  equivalent  electric  and  magnetic  currents  on  the  surfaces  of 
all  entities,  and  evaluates  the  net  electric  currents  in  all  probes. 
The  results  of  the  calculations  are  the  impedance  parameters 

(z)  [19]. 

To  model  the  measurements,  we  suppose  a  uniform  linear 
array  of  M  sensors  placed  in  front  of  the  reinforced  wall.  The 
array  takes  measurements  at  N  known  positions,  at  L  frequen¬ 
cies.  For  2D  modeling,  instead  of  induced  voltages  at  the  sensor 
ports  we  use  the  induced  electric  field  in  the  sensors.  The  mea¬ 
sured  electric  field  in  the  *th  sensor  is 
M 

Ei{n ,  fi)  =  Y  Eik(n,  fi )  +  u , 

k= 1 

i  =  1, . . . ,  M,  n  —  1, . . . ,  N,  l  =  1, . . . ,  L  (1) 

where  Eik(n,  fi )  is  the  induced  electric  field  in  the  ith  sensor 
when  the  /,:th  sensor  is  excited,  //  is  the  operating  frequency,  n 
is  the  index  of  the  measurement  position,  and  u  is  the  measure¬ 
ment  noise  and  modeling  error.  We  assume  complex,  zero-mean 
Gaussian  noise  with  variance  a2. 

The  antenna  array  can  be  considered  as  a  multiport  network. 
If  the  impedance  parameters  are  used,  the  measured  induced 
electric  field  reads 
M 

Ei(n,  fi)  =  Y  zik{n,  fi)Ik(fi)  +  u, 

k= 1 

i  =  1, . . . ,  M,  n  —  1 , .  , ,,,  N.  1  =  \ , L  (2) 

where  Zik(k,fi)  is  the  mutual  impedance  coefficient  between 
the  ith  and  fcth  sensors,  and  Ip,  is  the  feeding  current  of  the  /.  th 
sensor.  The  feeding  current  is  the  same  for  all  sensors.  We  as¬ 
sume  Gaussian  current  excitation  in  the  time  domain.  The  mag¬ 
nitude  of  the  pulse  is  1A,  and  it  is  centered  at  t  =  0. 


C.  Time-Domain  Response  of  Homogeneous  and  Reinforced 
Walls 

The  impulse  response  of  concrete  walls  with  rebar  is  inves¬ 
tigated  in  [14],  [20].  Until  the  signal  reflected  from  the  bars 
reaches  the  sensors,  the  reflected  signals  from  a  homogenous 
wall  (without  the  rebar)  and  from  the  reinforced  wall  with  the 
same  dielectric  permittivity  are  identical.  Due  to  the  multiple  re¬ 
flections  from  the  bars,  the  signal  reflected  from  the  reinforced 
wall  has  an  oscillating  nature.  These  oscillations  make  the  esti¬ 
mation  of  the  objects  behind  the  wall  more  difficult  compared 
with  the  case  of  a  homogeneous  wall.  The  reinforcement  does 
not  introduce  additional  delay.  We  illustrate  this  property  in 
Fig.  2,  where  we  show  the  signals  reflected  from  a  perfectly  con¬ 
ducting  plane  (PEC)  placed  behind  a  reinforced  wall  and  behind 
a  homogeneous  wall.  (The  reflection  from  the  PEC  is  obtained 
by  subtracting  the  response  of  the  wall  from  the  response  of  the 
wall  with  the  PEC  behind.) 

The  simulated  excitation  waveform  is  a  Gaussian  pulse 
g(t.)  =  e^~aty ,  where  a  =  3e-9  s-1.  The  pulse  duration  is 
1.2  ns  (pulse  duration  is  defined  as  the  time  interval  in  which 
the  signal  amplitude  is  at  least  5%  of  the  maximum  value). 
The  frequency  response  is  calculated  from  5  MHz  to  2  GHz 
in  5  MHz  steps.  (The  usage  of  very  low  frequencies  such  as 
5  MHz  is  not  necessary  for  the  through-the-wall  imaging  since 
the  reflection  coefficient  of  the  reinforced  wall  is  very  high  at 
those  frequencies.) 

We  assume  the  wall  thickness  is  w  =  0.2  m  and  the  wall 
length  is  8  m.  The  bar  period  is  d\Mr  =  0.15  m  and  the  bar 
diameter  is  D\mr  =  1  cm.  The  permittivity  of  wall  material  is 
£r  =  3  —  j0.15.  The  transmitting-receiving  pair  is  at  a  distance 
of  0.75  m  from  the  front  side  of  the  wall,  the  separation  between 
the  sensors  is  0.2  m,  and  the  distance  between  the  PEC  and  the 
sensors  is  1.5  m. 

III.  Estimation  of  the  Wall  Dielectric  Permittivity 
and  Thickness  Using  Beamforming 

A.  Wall-Permittivity  Estimation 

Electromagnetic  waves  are  reflected  from  interfaces  of  media 
with  different  electromagnetic  properties.  These  reflections  can 
be  separated  in  the  time  domain  if  the  sensors  transmit  suffi¬ 
ciently  short  pulses.  The  dielectric  permittivity  of  the  reinforced 
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Fig.  3. 
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Scheme  used  for  wall-permittivity  estimation. 
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Fig.  4.  Scheme  used  for  wall-thickness  estimation. 


wall  can  be  calculated  in  the  same  manner  as  the  permittivity  of 
the  homogeneous  wall  in  the  time  gate  that  ends  when  the  signal 
reflected  from  the  bars  reaches  the  sensors  [20] .  Examples  of 
permittivity  estimation  using  time-domain  techniques  or  radar 
are  given,  e.g.,  in  [11],  [12],  [21],  and  [22].  Detailed  analysis  of 
the  electromagnetic  properties  of  concrete  can  be  found  in  [23]. 

In  order  to  separate  the  reflections  from  the  front  side  of  the 
wall  from  the  reflections  of  the  rebar,  the  duration  (T)  of  the 
transmitted  pulse  should  satisfy  T  <  Tmax  =  2 d^fe^/c,  where 
d  is  the  distance  of  the  rebar  from  the  front  side  of  the  wall 
and  c  is  the  propagation  speed  of  the  electromagnetic  waves  in 
the  air.  In  the  considered  case,  d  =  w/2  and  Tmax  =  1.2  ns. 
(Alternatively,  the  expected  resolution  of  the  system  is  Ad  = 
0.5  cT/y^.) 

The  reflection  coefficient  of  the  TM  wave  at  the  air-dielectric 
interface  is 


m = (3) 

cos  8  +  v  er  —  sin2  8 

in  the  case  of  small  losses.  Without  loss  of  generality,  we  as¬ 
sume  that  the  angle  8  is  small;  hence  R(8)  ~  R( 0).  To  estimate 
the  permittivity,  we  coherently  add  signals  received  by  the  sen¬ 
sors  in  the  array.  The  time  delay  (jik)  between  the  v'th  and  A:th 
sensors,  according  to  Fig.  3,  is  given  by 

&ik  —  \/ Od  )  2  T  (t/w  W 

T  \J (.Xk  ^  T  (y w  Vk')~ 

T~ik  =  $ikl  Q)  (4) 

where  the  coordinates  (xw,yw)  give  the  position  of  the  reflec¬ 
tion  point  on  the  wall.  The  focused  waveform  reads 

!  N  M  M 

E0 (t)  =  ^  1  ^  1  ^  .  Eik(n,  t.  +  Tik)  (5) 

n= 1 2=1 k= 1 

where  En~(n,t)  is  the  time-domain  waveform  of  the  electric 
field  induced  in  the  fcth  sensor  when  the  vth  sensor  is  excited, 
and  n  is  the  index  of  the  measurement  position.  The  waveform  is 
shifted  in  the  time  domain  with  respect  to  Tlk  to  enable  coherent 
summation  (focusing)  of  the  received  signals.  Multiplication  by 
the  factor  \J  Sik{n)  ensures  that  all  received  signals  have  the 
same  attenuation,  corresponding  to  the  attenuation  at  distance 
s  =  1  111. 

We  calculate  the  estimate  of  the  reflection  coefficient  for 
small  8  as 


.R(O)  =  I  /  Eo(t)h(t)dt  ]  /  [  /  h2(t)dt 


(6) 


where  h(t  —  s/cq)  is  the  induced  electric  field  in  the  sensor  that 
is  separated  by  distance  s  =  1  m  from  the  transmitting  sensor 
(sensors  are  in  vacuum).  We  will  refer  to  /;(/.  )  as  the  reference 
pulse.  The  estimate  of  the  wall  permittivity  from  (3)  is 


v+m)  ‘ 


B.  Wall-Thickness  Estimation 

The  signal  reflected  from  the  back  side  of  the  wall  contains  in¬ 
formation  about  the  wall  thickness.  In  the  case  of  the  reinforced 
wall,  this  reflection  is  concealed  by  the  multiple  reflections  from 
the  bars.  To  uncover  it,  we  focus  the  measured  signals  in  the  time 
domain  with  respect  to  the  delay  of  those  pulses. 

The  signal  (Fig.  4)  transmitted  by  the  ith  sensor,  reflected 
from  the  back  side  of  the  wall,  and  received  by  the  fcth  sensor  is 
delayed  for  time 

Tik  (w)=  (( 2  yw  -  yi  -  yk )  /  cos  8 + 2wer  /  ^-(sin  8)2)  /cq 

(8) 

where 

8  =  arg  min  {  xk  -  as,--  ( 2  yw -yi-yk)t<m8- 
6  1 

2tt;sin  8/ y/er  —  (sin#)2!.  (9) 

The  summation  of  the  received  signals,  now  focused  with  re¬ 
spect  to  the  delay  Tik(w),  is 

.  N  M  M 

Eo(t)  =  Yl  Yl  Yl Eik  (n't+ rifcW)  •  (10) 

n= 1  2=1  k=  1 

The  wall  thickness  is  estimated  from  the  correlation  of  the 
focused  electric  field  with  the  reference  pulse  h(t) 

T 

I(w)  =  j  Eo(t)h(t)dt.  (11) 

o 

We  calculate  the  correlation  for  different  wall  thicknesses,  w. 
The  estimate  of  the  wall  thickness  is  defined  as  the  maximum 
of  (1 1).  Due  to  the  periodic  nature  of  the  reflected  field,  the  cor¬ 
relation  does  not  yield  a  unique  solution.  Nonetheless,  practical 
limits  on  the  wall  thickness  reduce  the  number  of  estimates  to 
only  a  few. 


IV.  Estimation  of  Target  Positions 

Movements  of  people  behind  the  wall  perturb  the  outside  an¬ 
tenna  measurements.  We  assume  the  measurements  of  the  sta¬ 
tionary  scene  behind  the  wall  are  available.  We  obtain  the  sig¬ 
nals  reflected  from  the  targets  by  subtracting  the  measurements 
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Fig.  5.  Position  of  sensors  and  target. 

of  the  stationary  background  from  the  measurements  altered  by 
the  appearance  of  the  people.  We  apply  beamforming  [  1]  on  the 
reflected  signals  to  estimate  the  positions  of  the  moving  targets. 
The  proposed  technique  can  be  applied  in  the  case  for  which  the 
characteristics  of  the  bars  are  unknown.  However,  we  show  that 
the  estimation  is  significantly  improved  if  we  have  some  prior 
knowledge  of  the  bars. 

A.  Estimation  When  Bar  Parameters  are  Unknown 

The  space  behind  the  wall  is  divided  into  small  cells  [1].  We 
focus  the  array  response  to  each  cell  to  determine  the  presence  of 
the  targets.  We  assume  the  speed  of  the  electromagnetic  waves 
in  the  reinforced  wall  to  be  the  same  as  the  speed  in  a  homoge¬ 
neous  wall  of  the  same  permittivity. 

If  the  target  is  located  in  the  cell  with  center  point  ( x ,  y),  the 
fcth  sensor  receives  the  signal  transmitted  from  the  ?'th  sensor 
delayed  for  Ty(x,  y)  +  rk(x ,  y),  where  n,  (ry.)  is  the  signal  prop¬ 
agation  time  from  the  ith  (7,:th)  sensor  to  the  target  (Fig.  5) 

Ti,k(x,  y)  =  {(y  -  y^k  -  w)/  cos  0;,fc+ 

wer/y4r  -  (sin^fc)2^  /cq.  (12) 

The  signal  transmitted  by  the  ith  sensor  is  launched  at  angle  8i, 
while  the  /.:th  sensor  receives  the  signal  reflected  from  the  target 
at  angle  6k.  Those  angles  are  calculated  as 

Oi,k  =  argmin  \x  -  Xi,k  -  (y  -  yi,k  ~  w)  tan  6- 

e  1 

w  sin  0/ \Je r  -  (sin#)2  j .  (13) 

We  coherently  add  the  received  signals  with  respect  to  the  delay 
Ti(x,  y)  +  Ty. (x.  y).  The  focused  electric  field  is 

N  M  M 

AE(t-,X,y)  =  £££  A Eik  ( n ,  t  +  Ti(x ,  y)  +  rk(x,  y)) 

n= 1  2=1  k= 1 

AEik  =Eik  -  Efkat  (14) 

where  i?|^at  is  the  induced  electric  field  in  the  case  of  the  sta¬ 
tionary  scene  and  Erk  is  the  induced  electric  field  due  to  the 
appearance  of  the  moving  objects.  The  image  pixel  at  (x,  y)  is 
obtained  by  correlating  the  focused  electric  field  with  the  refer¬ 
ence  pulse  h(i) 

T 

I(x,y)  =  J  AE(t;x,y)h(t)dt.  (15) 

o 

B.  Estimation  When  Bar  Parameters  are  Known 

The  reference  pulse,  h(t),  represents  the  waveform  of  the  in¬ 
duced  electric  field  in  the  receiving  sensor  for  a  given  excitation 


0  o 
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+  X  0  o 

xhy,  a 

A 

Fig.  6.  Position  of  the  sensors  for  numerical  calculation  of  the  wall  transmis¬ 
sion  coefficients. 


Fig.  7.  Transmission  coefficient  of  the  wall  as  a  function  of  the  sensor  positions. 
We  adopt  Xx  =  x2  =  0,  y i  =  —y2  =  —y  (Fig.  6). 

waveform.  Sensors  are  assumed  to  be  in  a  vacuum.  The  influ¬ 
ence  of  the  wall  on  the  received  waveform  is  not  taken  into  ac¬ 
count.  Here,  we  model  the  waveform  distortion  due  to  the  pres¬ 
ence  of  the  wall  with  bars.  We  assume  that  the  distance  between 
the  bars  and  their  diameter  are  known.  The  estimation  of  the 
unknown  bar  parameters  is  studied  in  [11],  [12].  The  authors 
developed  a  multiple  stratification  model  of  the  reinforced  con¬ 
crete,  using  the  characteristic  impedance  of  the  wire  net  derived 
by  Marcuvitz  [25]. 

The  transmission  coefficient  of  the  reinforced  wall  for  the 
plane  wave  incidence  is  given  by 

T(f,6)  =  Et/Ei  (16) 

where  0  is  the  incident  angle,  E\  is  the  incident  electric  field, 
and  Et  is  the  transmitted  electric  field.  There  are  several  studies 
of  the  transmission  and  reflection  coefficients  of  the  infinitely 
long  reinforced  wall  for  the  plane-wave  incidence,  e.g.,  [8],  [9]. 
For  the  problem  that  we  consider  here,  the  reinforced  wall  is  not 
necessarily  in  the  far  field  of  the  sensors.  Hence,  the  transmis¬ 
sion  through  the  wall  depends  on  the  position  of  the  sensors.  In 
order  to  compute  the  transmittivity  of  the  wall  in  the  near  field, 
we  use  the  model  shown  in  Fig.  6.  We  assume  that  the  sensor 
located  at  (x\ ,  y\ )  is  excited  and  compute  the  induced  electric 
field  in  the  sensor  located  at  ( X2,y 2).  We  define  the  transmis¬ 
sion  coefficient  as 

T(f,6)  =  E12(f,6)/E°2(f,6)  (17) 

where  E\  o  is  the  induced  electric  field  in  the  receiving  sensor 
with  the  wall  being  present,  and  E\}2  is  the  induced  electric  field 
when  the  sensors  are  in  a  vacuum.  (The  dependence  of  (17)  on 
X\ ,  7/1 ,  x2 ,  and  y2  is  assumed.) 

In  Fig.  7,  we  show  the  variation  of  the  wall  transmission  co¬ 
efficient  as  a  function  of  the  sensor  positions  for  two  different 
values  of  (4>ar/Ac  (Ac  is  the  wavelength  in  the  wall).  We  set  the 
sensor  positions  to  (xi  =  0,  y\  =  —y)  and  (x2  =  0,  y2  =  y). 
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Fig.  8.  Normalized  pulse  shapes  h(t)  calculated  for  various  incidence  angles. 
Reference  pulse  h(t )  is  denoted  by  line  without  markers. 


w[m] 


Fig.  9.  Estimation  of  wall  thickness:  correlation  of  focused  electric  field  and 
reference  pulse. 


In  the  vicinity  of  dbar/Ac  =  1,  the  transmission  coefficient  of 
the  reinforced  wall  is  periodically  dependent  on  the  distance  be¬ 
tween  the  sensors  and  the  wall.  This  result  is  in  agreement  with 
the  result  reported  in  [10],  where  the  periodicity  of  the  transmit¬ 
ting  characteristics  of  the  reinforced  wall  was  noted. 

The  distances  between  the  sensors  and  the  wall  affect  the 
amplitude  of  the  received  signal;  their  influence  on  the  wave¬ 
form  of  the  received  signal  is  negligible  (up  to  the  scaling  con¬ 
stant).  Since  we  are  interested  in  the  shape  of  the  received  wave¬ 
form,  we  compute  the  transmission  coefficients  at  fixed  sensor 
positions. 

The  signal  reflected  from  the  target  behind  the  wall  passes 
through  the  wall  twice.  For  the  target  from  Fig.  5,  the  signal 
transmitted  at  angle  6i  passes  through  the  wall,  gets  reflected 
from  the  target,  passes  again  through  the  wall,  and  reaches  the 
sensor  at  angle  6k-  The  induced  electric  field  in  the  receiving 
sensor  is  approximately 

Eik{f)~H{f)T{fA)T{f,6k)  (18) 

where  H(f)  is  the  Fourier  transform  of  Hence,  we  adopt 
for  a  new  reference  pulse: 

h(t)  =  rect  (^jF-\AH{f)T{f,  6k)  exp(j27r/r0)} 

(19) 

where  to  is  the  time  delay  that  assures  that  the  reference  pulse 
is  centered  at  t  =  0,  and  A  is  the  normalization  constant. 

In  Fig.  8  we  show  the  waveforms  induced  in  the  receiving 
sensor  after  reflecting  from  PEC  behind  the  wall  (inset  of  Fig.  2). 
The  waveforms  are  calculated  for  various  angles  of  incidence, 
normalized,  and  centered  at  t  =  0.  In  the  same  figure,  we  also 
show  the  reference  pulse  h(t),  denoted  by  the  solid  line.  Since 
the  pulse  shape  does  not  vary  notably  with  the  incident  angle, 
we  use  the  approximation  of  (19)  for  the  reference  pulse: 

h(t)  =  rect  F-1  j AH(f)  (T(/,0o))2exp(j27r/r0)}  . 

(20) 

V.  Results 

We  use  the  electromagnetic  solver  from  Section  II  to  simu¬ 
late  the  measurements  of  the  scene  behind  a  reinforced  wall. 
The  number  of  moving  objects,  as  well  as  their  properties,  is 
unknown. 


We  assume  a  wall  of  thickness  w  =  0.2  m  and  length  8  m. 
We  consider  various  bar  periods:  dbar  =  0.15  m,  dbar  =  0.1  m, 
and  ((bar  =  0.05  m.  The  bar  diameter  is  set  to  D\,Hr  =  1  cm. 
The  adopted  relative  permittivity  of  the  wall  material  is  £,-  = 
3  —  j0.15.  The  moving  objects  are  modeled  as  perfectly  con¬ 
ducting  cylinders  of  radius  r  =  0.2  m,  centered  at  x  =  1  m, 
y  =  1.2  m  and  x  =  —1.5  m,  y  =  0.75  m,  respectively.  Be¬ 
sides  the  targets,  there  are  static  (clutter)  objects  representing 
furniture  and  interior  walls.  We  modeled  the  static  objects  as 
dielectric  rectangles.  The  adopted  relative  permittivity  for  the 
static  objects  is  the  same  as  the  wall  permittivity. 

The  measuring  system  consists  of  a  uniform  linear  array  of 
M  =  3  sensors.  The  separation  between  adjacent  sensors  is 
0.4  m.  The  array  moves  parallel  to  the  wall  at  a  distance  of 
0.75  m  from  the  front  side  of  the  wall.  The  array  takes  measure¬ 
ments  every  0.2  m.  The  total  number  of  measurement  locations 
is  32. 

The  sensors  transmit  a  Gaussian  pulse,  g(t),  defined  in 
Section  II-C.  We  also  examine  the  cases  in  which  the  spectrum 
of  the  signal  is  shifted  in  the  frequency  domain,  i.e.,  G(f  —  /o), 
0  <  /o  <  1  GHz.  The  frequency  response  is  calculated  from 
5  MHz  to  2  GHz  in  5  MHz  steps. 

The  measurements  are  corrupted  with  white  noise.  We  cal¬ 
culate  the  signal-to-noise  ratio  with  respect  to  the  power  of  the 
electric  field  induced  in  the  sensor.  Most  of  the  induced  elec¬ 
tric  field  is  due  to  the  direct  coupling  between  the  sensors.  The 
power  of  the  electric  field  scattered  from  the  target  is  several  or¬ 
ders  of  magnitude  less  than  the  total  electric  field  induced  in  the 
sensor.  For  this  reason,  the  signal-to-noise  ratio  cannot  be  very 
low. 

A.  Unknown  Bar  Parameters 

In  this  example  we  set  dbar  =  15  cm  and  /q  =  0.  Using 
(3) — (7),  we  obtain  £r  =  3  for  the  estimate  of  the  wall  per¬ 
mittivity.  We  apply  (11)  to  compute  the  estimate  of  the  thick¬ 
ness.  The  result  is  shown  in  Fig.  9.  The  dotted  curve  represents 
the  correlation  of  the  reference  pulse  and  the  focused  electric 
field.  There  is  a  maximum  corresponding  to  the  real  wall  thick¬ 
ness,  w  =  0.2  m.  There  are  also  maxima  for  w  =  0.3  m  and 
w  =  0.4  m.  These  other  maxima  are  a  consequence  of  the  fo¬ 
cusing  reflections  from  the  bars.  The  separation  between  the 
maxima  depends  on  the  bar  period.  For  comparison,  we  repeat 
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Fig.  10.  Image  of  the  scene  behind  the  wall  computed  using  £r  =  3,  /o  =  0, 
and  (a)  w  =  0.2  m  (real  wall  thickness)  and  (b)  w  —  0.3  m  (erroneous  wall 
thickness).  The  waveform  distortion  is  not  taken  into  account.  The  rebar  period 
(dbar  =  15  cm)  is  assumed  to  be  unknown.  The  adopted  SNR  is  30  dB. 


the  analysis  for  the  case  in  which  the  bar  period  is  =  25  cm 
and  all  other  parameters  are  the  same.  The  result  is  also  shown 
in  Fig.  9  (solid  line).  Again,  we  have  a  maximum  corresponding 
to  the  real  wall  thickness,  w  =  0.2  m,  as  well  as  other  spurious 
maxima.  Hence,  the  beamforming  does  not  produce  a  unique 
value  for  the  wall  thickness. 

We  estimate  the  position  of  the  targets  using  the  wall  per¬ 
mittivity  and  thickness  estimates:  er  =  3,  w  =  0.2  m  and 
ir  =  3 ,w  =  0.3  m,  along  with  dbar  =  15  cm.  The  adopted 
signal-to-noise  ratio  is  SNR  =  30  dB.  (The  power  of  the 
field  reflected  from  the  targets  is  20  dB  less  than  the  power 
of  the  total  field  induced  in  the  sensors.)  The  image  calculated 
using  the  exact  wall  thickness  is  shown  in  Fig.  10(a),  and  the 
image  calculated  using  the  erroneous  wall  thickness  is  shown 
in  Fig.  10(b).  The  black  lines  denote  the  true  positions  of  the 
cylinders  and  clutter  objects.  Due  to  the  oscillating  nature  of 
the  reflected  signal,  the  images  are  blurred.  The  blurriness  is 
more  pronounced  when  the  erroneous  wall  thickness  is  used 
[Fig.  10(b)], 

B.  Modeling  Bar  Distortion 

When  the  bar  characteristics  are  unknown,  the  analysis 
reduces  to  the  case  in  which  the  objects  are  hidden  behind  a 
homogeneous  wall  of  the  same  permittivity.  However,  the  target 
spreads  calculated  in  this  way  may  be  large.  If  information 
about  the  bars  is  available,  we  can  refine  the  estimation  by 
modeling  the  influence  of  the  bars  on  the  signal  waveform.  We 
examine  this  improvement  for  various  bar  periods  and  different 
frequency  content  of  the  excitation  waveform. 

We  repeat  the  experiment  from  Fig.  10  assuming  the  bar  pa¬ 
rameters  are  known.  The  images  are  calculated  using  er  =  3, 
w  =  0.2  m  [Fig.  11(a)]  and  er  =  3 ,w  =  0.3  m  [Fig.  11(b)]. 
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Fig.  1 1.  Image  of  the  scene  behind  the  wall  computed  using  er  =  3,  f0  = 
0,  and  (a)  w  =  0.2  m  (real  wall  thickness)  and  (b)  w  =  0.3  m  (erroneous 
wall  thickness).  The  waveform  distortion  is  taken  into  account.  The  rebar  period 
(dbal.  =  15  cm)  is  assumed  to  be  known.  The  adopted  SNR  is  30  dB. 


The  image  quality  [Fig.  11(a)]  is  significantly  improved  when 
the  corrected  pulse  shape  h(t)  is  used  compared  with  the  case  in 
which  the  distortion  is  not  modeled  [Fig.  10(a)].  The  pixels  with 
the  most  intensive  colors  now  clearly  represent  the  contours  of 
the  objects  (as  seen  by  the  sensors).  The  image  obtained  using 
erroneous  wall  thickness  is  again  more  blurred  [Fig.  11(b)], 

We  further  study  the  influence  of  the  SNR  on  the  estima¬ 
tion  accuracy.  We  compute  the  image  of  the  same  scene  for 
SNR  =  10  dB,  which  is  anticipated  as  the  signal-to-noise 
threshold  in  this  example.  The  target  traces  are  clearly  visible  if 
the  bar  parameters  are  known  [Fig.  12(a)].  In  contrast,  for  low 
SNR,  the  targets  cannot  be  discerned  without  waveform  cor¬ 
rection  [Fig.  12(b)]. 

Transmission  and  reflection  coefficients  of  the  reinforced 
wall  are  greatly  influenced  by  the  ratio  of  the  signal  wavelength 
and  the  bar  spacing,  [8],  [9],  [11],  [12].  The  investigations  in 
[11],  [12]  showed  that  when  the  ratio  dbar/Ac  is  decreased  to 
0.3-0. 4,  the  transmittivity  of  the  wall  is  significantly  reduced 
due  to  the  rebar.  Conversely,  the  influence  of  the  bars  is  reduced 
if  dbar/Ac  is  larger  than  0.7-0. 8.  We  study  those  results  in  the 
context  of  the  through-the-wall  imaging. 

As  the  central  frequency  of  the  pulse  spectrum  increases,  the 
position  of  the  targets  can  be  estimated  even  in  the  case  for 
which  the  bar  parameters  are  unknown.  Our  simulations  showed 
that  in  the  considered  example  (d\mr  =  15  cm,  SNR  =  10  dB), 
without  prior  knowledge  of  the  bar  parameters,  the  minimal  cen¬ 
tral  frequency  that  allows  estimating  targets  is  /(l  =  0.5  GHz. 
This  finding  is  in  agreement  with  the  results  from  [11],  [12], 
that  showed  that  the  attenuation  is  significant  for  all  frequen¬ 
cies  smaller  than  0.5  GHz  for  given  wall  parameters.  In  Fig.  13 
we  show  the  results  obtained  without  the  correction,  for  /o  = 
0.5  GHz  and  f0  =  1  GHz. 
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Fig.  12.  Image  of  the  scene  behind  the  wall  computed  for  SNR  =  10  dB, 
£r  =  3,  w  =  0.2  m,  dbar  =  15  cm,  and  f0  =  0.  The  waveform  distortion  is 
(a)  taken  and  (b)  not  taken  into  account. 


Fig.  14.  Image  of  the  scene  behind  the  wall  computed  for  SNR  =  15  dB, 
£r  =  3,  w  =  0.2  m,  dbar  =  10  cm,  and  /o  =  0.75  GHz.  The  waveform 
distortion  is  (a)  taken  and  (b)  not  taken  into  account. 


(b) 


Fig.  13.  Image  of  the  scene  behind  the  wall  computed  for  SNR  =  10  dB, 
£r  =  3,  w  =  0.2  m,  dbar  =  15  cm,  and  (a)  /o  =  0.5  GHz  and  (b)  and 
/0  =  1  GHz.  The  waveform  distortion  is  not  taken  into  account. 


Fig.  15.  Image  of  the  scene  behind  the  wall  computed  for  SNR  =  20  dB, 
£r  =  3,  w  =  0.2  m,  dbar  =  0  cm,  and  /0  =  1  GHz.  The  waveform  distortion 
is  (a)  taken  and  (b)  not  taken  into  account. 


We  repeated  the  experiment  for  dbar  =  10  cm  and  SNR  = 
15  dB  (Fig.  14).  The  improvement  due  to  the  waveform  cor¬ 
rection  is  significant  for  /o  <  0.75  GHz,  while  there  is  almost 
no  difference  at  /o  =  1  GHz.  According  to  [11],  [12],  in  this 
case,  the  rebar  has  major  influence  for  frequencies  smaller  than 
0.7  GHz. 

Finally,  for  ribar  =  6  cm,  it  was  expected  that  the  im¬ 
provement  would  be  significant  for  /o  <  1  GHz.  (For 

dbar  =  6  cm  the  attenuation  is  considerable  for  frequencies 
smaller  than  1  GHz.)  Our  computations  confirm  this  conclu¬ 


sion.  In  Fig.  15  we  show  the  images  computed  for  fo  =  l  GHz 
and  SNR  =  20  dB,  with  and  without  correction. 

We  examined  the  robustness  of  the  algorithm  to  the  errors 
in  bar  parameters.  In  Fig.  16  we  show  the  results  obtained  as¬ 
suming  that  the  bar  period  was  ribar  =  10  cm,  while  the  true 
value  was  <7bar  =  15  cm.  The  image  does  not  differ  significantly 
from  the  image  obtained  with  real  wall  parameters  [Fig.  12(a)], 

VI.  Conclusion 

We  addressed  the  important  problem  in  urban  warfare  of  es¬ 
timating  moving  targets,  such  as  personnel,  behind  a  reinforced 
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Fig.  16.  Image  of  the  scene  behind  the  wall  computed  for  SNR  =  10  dB, 
e,.  =  3.  in  =  0.2  m,  dbai-  =  11  cm,  and  f0  =  0  GHz.  Waveform  distortion  is 
modeled  using  erroneous  bar  period  (dbal.  =  10  cm). 


wall,  using  radar  measurements.  Reinforced  walls  significantly 
attenuate  low  frequencies  and  distort  the  transmitted  wave¬ 
forms.  The  signals  reflected  from  the  objects  behind  reinforced 
walls  are  oscillating  in  nature  and  are  of  long  duration  because 
of  the  bar  periodic  arrangement.  First,  we  considered  the  case 
in  which  the  reinforced  wall  is  completely  unknown.  We  added 
coherently  measured  electric  fields  in  various  time  gates  to 
estimate  wall  thickness  and  concrete  permittivity  indepen¬ 
dently.  We  assumed  the  availability  of  the  measurements  of 
the  static  scene  behind  the  wall  (walls,  furniture,  etc.)  and 
measurements  altered  by  the  appearance  of  people.  The  differ¬ 
ence  between  these  two  measurements  is  associated  with  the 
electromagnetic  field  scattered  from  the  moving  targets.  We 
applied  beamforming  to  the  field  reflected  from  the  targets  to 
estimate  their  number  and  locations.  The  algorithm  performed 
satisfactorily  in  the  case  in  which  the  bar  diameter  and  period 
are  both  unknown.  Flowever,  due  to  multiple  echoes,  the  target 
spreads  are  large.  We  also  analyzed  the  problem  in  which  the 
characteristics  of  the  metallic  rebar  are  known.  We  improved 
the  estimation  significantly  by  modeling  the  waveform  dis¬ 
tortion  due  to  the  bars.  The  resulting  images  are  focused  and 
accurately  represent  the  contours  of  the  targets.  The  algorithm 
is  robust  to  the  ambiguities  in  bar  parameter  values.  In  addition, 
the  minimal  necessary  SNR  is  lower  compared  with  the  case  in 
which  the  influence  of  the  bars  on  the  signal  shape  is  ignored. 
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